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Abstract 

6 

^ We study a cosmological scenario in which the DBI action governing the motion of 

a D3-brane in a higher-dimensional spacetime is supplemented with an induced gravity 
term. The latter reduces to the quartic Galileon Lagrangian when the motion of the 
^ brane is non-relativistic and we show that it tends to violate the null energy condition 

and to render cosmological fluctuations ghosts. There nonetheless exists an interesting 
^ parameter space in which a stable phase of quasi-exponential expansion can be achieved 

O while the induced gravity leaves non trivial imprints. We derive the exact second-order 

00 action governing the dynamics of linear perturbations and we show that it can be simply 

2 understood through a bimctric perspective. In the relativistic regime, we also calculate 

^ the dominant contribution to the primordial bispectrum and demonstrate that large non- 

>■ Gaussianities of orthogonal shape can be generated, for the first time in a concrete model. 

^ More generally, we find that the sign and the shape of the bispectrum offer powerful 

^ diagnostics of the precise strength of the induced gravity. 
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1. Introduction 

On-going measurements of the cosmic microwave background (CMB) anisotropies by the 
Planck satelhte [1] promise to provide us with yet more precise information on the very 
early universe than those we gained from the WMAP satellite [2]. In this respect, while 
Gaussianity of the primordial fluctuations is up until now still preferred by the data [3, 
4], a tremendous amount of attention is currently paid to their possible deviation from 
perfect Gaussian statistics (see [5, 6] for recent observational reviews). Indeed, any such 
detection would shed some light on the interactions of the field (s) driving inflation or 
its alternatives, a piece of information that can not be probed by the Gaussian linear 
theory. These interactions are very small in slow-roll single field inflationary models, which 
therefore predict a very small amplitude of primordial non-Gaussianities [7]. However, 
going beyond these simplest scenarios have demonstrated that sizeable non-Gaussianities 
can be generated in a wide variety of early universe models (see for instance the reviews 
[8-11] and references there in). 

Amongst them, the Dirac-Born-Infeld (DBI) model [12-16], in which the inflaton is 
identified with the radial coordinate of a probe D3-brane in a warped geometry, has at- 
tracted considerable amount of attention as an inflationary scenario inspired by string the- 
ory, producing large non-Gaussianities, and possibly alleviating the famous eta-problem 
[17, 18]. It is also one of the few /c-inflationary models [19, 20] whose higher derivative 
structure is known not to be spoiled by quantum corrections [21]. One of the lessons 
we learnt from the extensive studies to put constraints on DBI inflation [22-32] is that 
so-called ultraviolet models are under strain from observations because they generate a 
too high level of non-Gaussianities^. However, this applies only to single field models of 
DBI inflation, in which the angular degrees of freedom of the brane are assumed to be 
frozen. As was shown in Ref. [34] and in further details in [35], relaxing this assumption 
and allowing the brane to evolve and fluctuate in the angular directions (see e.g [36-40]) 
alleviates these constraints. Indeed, in such a multiple field scenario, the amplitude of 
equilateral non-Gaussianities is suppressed as a result of the transfer from entropic to cur- 
vature perturbations (see [41-50] for subsequent studies on multi-field inflationary models 
with non-canonical kinetic terms). 

Recently, de Rham and Tolley proposed a very interesting extension of the DBI action 
[51] which generalizes and unifies single field DBI and the modified gravity model called the 
Galileon [52-54] (see the work by Horndeski [55] for an early study of a related more general 
class of models and Refs. [56-76] for some studies of cosmology based on the Galileon field). 
Following the geometrical picture of a D3-brane embedded in a five-dimensional spacetime 
(see also [77-80]), they showed that only a finite number of terms are allowed that give 



Concerns also exist about their consistent embedding in string theory [25, 31, 33]. 
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rise to second-order equations of motion for the brane position modulus, and that these 
terms reduce to the so-called Galileon Lagrangians in the non-relativistic limit in which the 
motion of the brane is slow (see [81] for a study of the non-Gaussianities generated by one 
of these terms). This relativistic extension of the Galileon model has later been generalized 
to the case of an arbitrary number of extra dimensions in Ref . [82] , based on the work [83] 
that studied matching conditions for distributional sources of arbitrary co-dimension in the 
context of Lovelock gravity (see Refs. [84-95] for other interesting generalizations in the 
single field case and their cosmological implications). Multifield extensions of the Galileon 
terms were extensively studied in Refs. [96-101] but the cosmology of the multifield DBI 
Galileon model started to be explored by one of us only very recently [102]. In this paper, 
we give further details on this model formulated in an ambient ten-dimensional spacetime. 
In this CclSG clS relevant to string theory - only an induced gravity term can be added to 
the DBI action [82, 83]. As we will see, this geometrical structure implies that it is useful 
to think of this model not as a modification of multifield DBI inflation, but rather as a 
multi scalar-tensor modification of gravity with a peculiar bimetric structure. 

The layout of this paper is the following. In section 2 we introduce our model and 
derive the conditions under which it can sustain a phase of quasi de-Sitter expansion. 
We study the scalar and tensor linear perturbations about a cosmological background in 
section 3 while section 4 is devoted to the study of the non-Gaussianities of the primordial 
curvature perturbation. We give our conclusions in section 5 and leave some technical 
details to several appendices. 

2. Set-up and homogneous evolution 

2.1. Action and equations of motion 

The geometrical set-up we consider is the same as in standard multifield DBI inflation [35]: 
we consider a D3-brane with tension T3 evolving in a 10-dimensional geometry described 
by the metric 

ds^ = h~^/\y^) gf.^dx^'dx" + h^^\y^) Gu{y^) dy^dy^ = HAsdY^dY^ (2.1) 

with coordinates Y"^ = [x^,y^}, where = 0, . . . 3 and / = 1, . . . , 6 (the label / has been 
chosen in this way as below it will label the multiple effective scalar fields). An important 
role will be played by the induced metric on the 3-brane, 

= HABd.Y^i^d.Y^^) (2.2) 

where the brane embedding is defined by the functions Y'^^^(x^), with the x^ being the 
spacetime coordinates on the brane. In our case, they coincide with the first four bulk 
coordinates, so that, on writing Y|^^ = (x'^, (y9^(a;^)), we find 

7/.. = h-^'^ {g^, + h Gijd^^'d^) . (2.3) 
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As discussed in the introduction, our starting point four-dimensional effective action will 
be the sum of three terms: 

• the DBI action, whose kinetic part, proportional to the world- volume encompassed 
by the brane^, is simply a cosmological constant in terms of the induced metric: 

• the Einstein-Hilbert Lagrangian associated to the induced metric: y/—'jR['j]. This 
term is absent in standard DBI inflation and will be the focus of our attention^. 

• the standard Einstein-Hilbert Lagrangian associated to the cosmological metric: 



The total action thus reads 
S= I d^x 



2 V ^ L^j 2 
where Mp and M are constant mass scales^ and 



Ml 



(2.4) 



/:brane = " (v^ " l) " V{<P') . (2.5) 

Here, following the literature on brane inflation, we have introduced the rescaled variables 

f = ^. ¥ = V^/ , (2.6) 

we have included potential terms in the brane action and we have defined 

V = det(5^: + / Gu9^'dpct>'d,ct>') (2.7) 

where Gij{(f)^) will play the role of a metric in the space of the scalar fields (p^ . On defining 
the mixed kinetic terms for the scalar fields 

X''^-^^g^'^d,<p'd,<p', (2.8) 

^ Strictly speaking, the DBI action also depends on bulk fields as well as on the gauge field living on the 
brane. Their cosmological consequences were analysed in Ref. [43]. 

^ Such an induced gravity term is present in bosonic string theory [103]. It is absent at tree level in type 
II superstring theories [104] but may be induced at 1-loop string level. We thank Liam McAllister for 
useful discussions on this point. 

^ If some matter was to be added to the action, we have in mind that it would be minimally coupled to the 
metric g^^. Hence, we conventionally refer to the mass scale associated to the standard Einstein-Hilbert 
Lagrangian as Mp. However, we should bear in mind that, because of the induced gravity, this is a 
priori not directly related to the Newton constant as would be measured in laboratory experiments. 
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it has been shown in [35] that the exphcit expression of V reads 

V = l- 2fGijX^-^ + AfxfXj^ - SfxfxjX^^ + 16fxfxjX^xi\ (2.9) 

where the brackets denote antisymmetrisation on the field indices and X/ = GikX^^ . 
Similarly, one can express a/— 7i?[7] in terms of the fields and the geometrical quantities 
associated to the cosmological metric, leading to a multifield relativistic extension of the 
quartic Galileon Lagrangian in curved spacetime. The resulting expression, that can 
be found in appendix (A), is very intricate and rather obscures the physics which is 
at play by breaking the symmetry between the cosmological and the induced metric in 
the "gravitational" part of the action (2.4). In the following, we rather try to keep this 
symmetry manifest by treating the induced gravity and Einstein-Hilbert action on equal 
footing. 

Following this bimetric perspective, one can write the gravitational equations of motion 
- obtained by varying the action (2.4) with respect to g^j_y - in the compact form 



MlG^^lg] + M^I^G^^H = T,r_ (2.10) 
where we have used the fact that 

^ = h-^^g^, (2.11) 

we have defined 

T/^- = 2 ^(V^^brane) ^ _ ^ ^^M- _g^^^(v--\ (2.12) 

and ^^'^ denotes the inverse of the metric 

9f.. = 9>.. + fGjjd^(l)'d,(l)' (2.13) 

(the latter is simply h^^'^'j/j,^). As for the equations of motion for the scalar fields, they 
read 

+ j^^) {S^j + 2fGj,X-<) (n^. + t ) " S ("^^ + /^) = ° • 

(2.14) 

Here, X'^-^ = -lg''^d^(j)^df,(j)'^ and 

^ + f i^V^^^V,^^ (2.15) 

where denotes the covariant derivative associated to g^u and T^ab the Christoffel 
symbols associated to Gjj = fGij. Note that, because of the induced gravity, second- 
order derivatives both of the metric and of the scalar fields enter into the gravitational 
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(2.10) and scalar fields (2.14) equations of motion, a property whose consequences have 
recently been explored in a slightly different context [84, 85]. 



It turns out that the warping - the fact that h (and hence /) depends non-trivially on 
the scalar fields - complicates the understanding of the role of the induced gravity in a 
non essential way. Hence, in the remaining of the main body of this paper, we consider a 
constant warp factor, in which case our starting point action simplifies to 



S 



/=cst 



d^X 



Ml 



-9R[9] + 



(2.16) 



where = M"^ jyh is constant. For simplicity, we also refer to the metric (2.13) as the 
induced metric. The consequences of a non-trivial warping are addressed in appendix (E). 



2.2. Background evolution 

In this section, we specialize the above equations to homogeneous configurations of the 
fields in a flat Friedmann-Lemaitre-Robertson- Walker (FLRW) background spacetime, of 
metric 

g^.^dx^'dx" = -dt^ + a^{t)dx^ . (2.17) 
Like in standard brane inflation, it is useful to introduce the background value of V (2.9): 

cl = l-f&\ (2.18) 

where a = \J Gij(j)^(f)-^ plays the role of an effective collective velocity of the flelds^. The 
value of c|, distinguishes the slow-roll regime where c|, ~ 1 and the brane action (2.5) 
takes a canonical form, from the so-called relativistic, or DBI, regime, where the brane 
almost saturates its speed limit and c|, ^ 1. In that case, the whole non-linear structure 
of the action (2.5) must be taken into account and large non-Gaussianities are expected. 
It is in this phenomenologically more interesting regime that we are particularly interested 
in. 

To derive the background equations of motion, it is not necessary to use the explicit 
expressions in appendix (A) for the geometrical quantities associated to the induced metric, 
like R[g\. It is rather more useful to notice that, in the background, the induced metric 
(2.13) takes the form of a flat FLRW metric whose scale factor is the cosmological one, a, 
but whose cosmic time t is such that 

di=cvdt. (2.19) 



^ In this paper, we use the symbol c|,, and not the more conventional c^, because, contrary to standard 
brane inflation, this does not coincide with the speed of sound for (scalar) perturbations. 
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Introducing the corresponding Hubble and "deceleration" parameters, given respectively 
by ^ = -§ = — and 

a at c-p 

I dH , , 

e = -^ ^ = e + s, 2.20 

where e = —H/H^ and s = cv/{Hcx)), it is then straightforward to derive the expressions 
of geometrical quantities associated to the induced metric from well known results. For 
example, from Ru[g] = —3H^{1 — e), we deduce that Ru[g] = —3/7^(1 — e), and hence, 
using the transformation law of tensors: Rtt[g] = — 3c|,iJ^(l — e) = — 3if^(l — e — s). 
Following this method, the modified Friedmann equations are readily found to be 

,9 1/1 



3i72(M|+-^j =pbranc = ^+y -1) (2.21) 



, , M^H^ ( 3 / 1 \\ ^2 , ^ 



We do not reproduce here the background equations of motion for the fields because they 
take a complicated and not very illuminating form. They can be found in full generality 
in appendix (E). As noticed above, second-order derivatives of the fields enter into the 
gravitational equation of motions (and vice- versa). However, solving the associated system 
of linear equations is not required to understand interesting properties of the homogeneous 
evolution. In particular, the way we derived Eq. (2.22) makes it clear that second-order 
derivatives of the fields enter into it only through the parameter s. Hence, in a quasi 
de-Sitter inflationary spacetime in which the time evolution of every quantity is slow with 
respect to that of the scale factor - i.e. X / (HX) <^ 1 - s is not expected to qualitatively 
alter the solutions of (2.22) and thus can be neglected. Note also that the combinations 
of mass scales that enter into Eqs. (2.21) and (2.22) are respectively Mp + M'^/c^ and 
Ml + W^/cv Hence, because of the non-zero velocity of the fields - cj) ^ 1 - the effect 
of the bimetric structure is not a simple rescaling of Mp. A parameter that will prove to 
be useful to measure this effect is 

K = — (2.23) 

that is comprised between c|, and 1. Combining Eqs. (2.21) and (2.22), we find (neglecting 

/ - 1 \ (2.24) 



2k {cvfV +l-cv) 2 {cvMl + m) Vl 

From the second term, one deduces that, for generic values of M, e acquires a large negative 
contribution when c|, -C 1 , which tends to violate the null energy condition and to hamper 
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inflation. Barring cancellations between the positive and negative contributions to e in 
(2.24), one thus finds that achieving a phase of quasi de-Sitter expansion in the relativistic 
regime - c|, ^ 1 - requires, in addition to the usual condition to enter into the DBI regime 

cvfV > 1 , (2.25) 

that 

< 4m|. (2.26) 

(and hence k ~ 1). At this point, given this restrictive condition on the mass scale 
associated to the induced gravity, it seems hopeless to expect anything new from multifield 
DBI Galileon infiation compared to standard multifield DBI infiation. However, this guess 
is not correct. To see this, let us introduce the dimensionless quantity 

a.i^. (2.27) 

Combining this definition with Eq. (2.21), one finds that, in our regime, 3a ~ jM-^{cvfV) 
and hence, being the product of a small parameter by a large one, is not necessarily small. 
The approximate equation deduced from (2.24): 

;i-3a) (2.28) 



then shows that the induced gravity can have a non negligible effect on the background 
evolution. Note also that this contribution to e is always positive in the regime where the 
theory is ghost-free, as will be demonstrated by the analysis of cosmological perturbations 
in the next section. 



3. Linear cosmological perturbations 

In this section, we study the dynamics of linear perturbations about a homogenous cosmo- 
logical solution. For this purpose, one can envisage several approaches. One is to perturb 
up to first order the equations of motion. This is not suitable for quantizing the fiuctua- 
tions but this is rather straightforward and one easily deduces from this the fact the two 
Bardeen potentials differ. In other words, there exists an effective anisotropic stress in the 
multifield DBI Galileon model, whose consequences for the growth of perturbations and 
in the late universe would be interesting to study. This is shown in appendix (B). The 
modern approach to cosmological perturbations during infiation, pioneered by Mukhanov 
and Chibisov [105] and Sasaki [106] and further developed by Maldacena [107], is rather to 
use the ADM formalism [108] to calculate efficiently the action at second order in the per- 
turbations. From this action, one can easily derive the equations of motion for the linear 
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perturbations but also determine the normalization of the vacuum quantum fluctuations 
and hence the amplitude of primordial cosmological perturbations. A priori, one can ap- 
ply this strategy to the explicit form of the action given in Eq. (A. 14). However, this 
expression is complicated and the resulting calculations are very cumbersome. Similarly 
to our derivation of the background equations of motion in section 2.2, the more efficient 
technique we will employ is to treat, within the ADM formalism, the cosmological and the 
induced metric on equal footing. 

Let us finally note that, although we will be ultimately interested in quasi de-Sitter 
backgrounds in the relativistic regime, our following analysis is exact and holds irrespective 
of the details of the homogenous dynamics. 

3.1. Second-order action for the perturbations 

We use the ADM formalism in which the metric is written in the form 

g^^dx^dx'' = -N'^dt^ + hij {dx' + N'dt) {dx^ + N^dt) (3.1) 

where hij is the spatial metric, A^* is the shift vector and A^ the lapse function. As discussed 
in the introduction of this section, following our bimetric perspective, it proves to be useful 
to introduce their induced gravity counterparts, such that 

g^^dx^dx" = -N'^dt^ + hij{dx' + N'dt){dx^ + N^dt) . (3.2) 

From Eqs. (2.13) and (3.1), their explicit expressions are given by 

h^j = hi, + fGudi(t)'d,<p^ (3.3) 
AT, = h,,N^ = K^N^ + fGiA'd.ct)' = Ni + fGjA'd.ct)' (3.4) 
= - fGjj<p^<j)'^ + h'^NiNj - h'^NiNj (3.5) 

where h^^ is the inverse of hij and (p^ is fully inhomogeneous. The "gravitational" part of 
the action (the first two terms in Eq. (2.16)) can then be simply expressed as two copies 
of the Einstein-Hilbert action in the ADM form^: 

^grav = ^ / dtd^xv^ (^Ari?(3) + j^{E,,E'^ - E'')^ 

+ ^ J dt&'xVl (^NR^^^ + ^{E^jE'^ - E^ij (3.6) 

^ When going from Eq. (2.16) to Eq. (3.6), we discarded some boundary terms, corresponding to the 
Gibbons-Hawking- York boundary action [109, 110] and its induced-gravity counterpart. The relevance 
of this kind of terms for Gahleon-type interactions was recently discussed in Refs. [Ill, 112] in a shghtly 
different context. 
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where h = det{hij), R^^^ is the scalar Ricci curvature built from hij, and the symmetric 
tensor Eij is defined by 

Eij = \K - ^V.iV, - U/,N, (3.7) 

where V represents the covariant derivative associated to the metric hij (similar definitions 
hold of course for quantities without"). It is then straightforward to derive the constraint 
equations for the lapse and the shift. Because the specific form of the brane action in (2.16) 
is not important in this calculation, and because it actually simplifies the presentation, 
we derive these constraint equations for a general scalar field Lagrangian of the form 
P(X^"^, </)^), where the kinetic terms X^'^ , defined in (2.8), take the ADM form 

X'' = ^v'v' - h^W,<p'd,<f>' (3.8) 

with 

= - N'di<p^ . (3.9) 
The variation of the action with respect to then yields the energy constraint 

ill / ~ ~ ~ ~ \ 

2{N'^P-P^j)V^V^)+Ml {N^R^^'^ +E^- E,jE'^)+M^d i^N^R^^'^ + E^ - EijE'^j = 

(3.10) 

while the variation of the action with respect to the shift Ni gives the momentum constraint 



+ 2 



(3.11) 



where P{ij) = | {dP/ dX^"^ + dP/ dX'^^) . The remaining task is as follows: we decompose 
the lapse and shift as 

N = 1 + SN, Ni = f3^i + Vi (3.12) 

where d^Vi = 0; we solve the linearized constraints equations for the auxiliary variables 6N, 
f3 and Vi, insert their solutions back in the action and deduce the quadratic action in terms 
of the true propagating degrees of freedom: two tensor modes and scalar perturbations 
(there is no active source of vector modes). At linear order, these two types of fiuctuations 
are decoupled and can be treated separately. We consider scalar perturbations in the 
next two subsections and then move on to discuss gravitational waves. In practice, the 
intermediate calculations for the scalar sector can be simplified by making appropriate 
choices of gauge. In this respect, before discussing the multifield situation, it is instructive 
to consider the case of a single inflaton field (A^ = 1, Gu = 1), to which we now turn. 
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3.1.1. Single field case 

A particularly convenient gauge in the single field case is the uniform inflaton gauge in 
which the (single) scalar perturbation C, appears in the spatial metric hij in the form 
hij = aitYe^^dij while the inflaton is homogeneous = (pit). The induced quantities (3.3), 
(3.4), (3.5) then considerably simplify to hij = hij, Ni = Ni and N'^ = N"^ — 1 + c|,. 
Details of the calculations for a general ^-inflationary Lagrangian of the form P(X, 0), 
with X = —^g'^'^d^(f)du(f), can be found in appendix (B). The resulting second-order (scalar) 
action is of the form 



5(2) = / dtd^xa^ ( A{t)C^ - S(t)^ ) , (3.13) 



a2 



which simply expresses the fact that ( is exactly massless and hence conserved on large 
scales - in case the decaying mode is indeed decaying - (see [113, 114] for a non-linear 
proof in a broader context). Concentrating on the DBI Galileon case, A and B can be 
cast in the following simple form: 



Ait) = ^{eK' + 34(1 - k"")) + ^ (iK^ + 34 fl - ^"j V (3.14) 

Bit) = M|(eK(3K-2) + 1) + ^ (^iK (^3-^ -2^ + K- cl^ , (3.15) 

where we remind the reader that e, defined in Eq. (2.20), is the deceleration parameter 
associated to the background induced metric and k, was defined in Eq. (2.23). Hence, 
this formulation reveals how the fluctuation ( reacts to the two background geometries: 
without the induced gravity, = 0, k = 1 and one recovers the standard fc-infiationary 
result [20] with a speed of sound Bit) /Ait) = 4- ^^e other extreme limit, if one 
formally considers Mp = 0, then k = 4 ^ind one finds a similar expression 

^(2),M-o = / did^xa^iM^cl (^(ly - I (3.16) 

in terms of the quantities associated to the induced metric with the replacement 4 ~^ ^v^- 
This can be easily understood: one can then treat the metric g as the cosmological metric^, 
in which case the action for gravity becomes canonical while the brane action can be 
expressed as Sbrane = / d'^X\/^-gPiXg,(l)) where Xg = —\g^''d^(t)dy(l) and 

P=(^--V^ ^l + 2fXg~j. (3.17) 



^ This argument is not completely straightforward and can not be applied blindly to the multifield case 
for example. Details about its validity can be found in appendix (D). 
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Using the general expression derived in [20], it is easy to verify that the speed of sound 
is given by l/c|, for the Lagrangian (3.17), with the consequence that the speed of sound 
with respect to the cosmological cosmic time is unity. In the general situation, one can not 
define an Einstein frame in which the gravitational action becomes canonical and one must 
resort to the full calculation, leading to the result (3.13) that nicely interpolates between 
the two extreme cases aforementioned, the relative importance of the two geometries being 
determined by the parameter n. However, the form (3.14) of the kinetic term is not 
appropriate for discussing the possible presence of a ghost. Because the conditions for 
avoiding the presence of ghosts are, a priori, more restrictive in the multifield case, we 
now turn to this situation, restricting our attention to two fields (/ = 1,2) and a trivial 
field space metric Gjj = 5ij for simplicity of presentation. 

3.1.2. Multifield case 

Contrary to the previous subsection, the fluctuations of all the scalar fields can not be 
erased in the multifield case by a suitable choice of time-slicing. Instead, a convenient 
gauge is the spatially flat gauge, in which the spatial metric takes its unperturbed value - 
hij = a'^{t)Sij - while the scalar degrees of freedom are directly taken to be the scalar fields 
fluctuations such that 0^ = (f)^ {t)+Q\t, x) (we omit the bar on (p^ in the following when 
there is no ambiguity). Details of the linearized constraint equations and their solutions 
can be found in appendix (B). Intermediate steps in the calculation of the second-order 
action are very lengthy and not illuminating. However, the final exact result can be cast 
in a very useful and elegant form. For that purpose, it is convenient to decompose the 
scalar field fluctuations into 



the background trajectory in field space and is the unit vector orthogonal to e^. In this 
way, the so-called adiabatic perturbation Q^r inherits all the properties of the singe field 
case while Qs, called the (instantaneous) entropy perturbation, embodies the genuinely 
multifield effects [115] (there would be — 1 entropy perturbations in general). Indeed, 
entropy fluctuations are fluctuations away from the background trajectory, something that 
can not happen in the single field case in which fluctuations only represent time-delays 
along the background trajectory. 

After going to conformal time r = J dt/a{t), and, upon using the canonically normal- 
ized fields = and Vg = w Qs with 



where = 0^ 




a a 



((1 — 9a)/t^ -|- 6aKc' 



.2 N 1/2 



(3.19) 



z = 




w = 




1/2 



(3.20) 




14 



the second-order action can be expressed in the remarkably simple form: 

■5(2) = lf dTcl'iju;^ + - 2^K>, - cUdv,)^ - c^.idv,)" - 2cl,dv,dv, 



+7"- + 



w" 
w 



(3.21) 



where a prime denotes a derivative with respect to conformal time and 



cl = ((1 -9a)K + 6a4) ^ 
1 + 



K — C 



((1 — 5q;)k + 2a;c|,) + Qasc^- — 



,2 n 



r2 



K 



2sa K — 
1 — 3a 1 — c|, 
c|,(l-«:)(l-c|,) 



2Mj,HVs 



a3 



(/t(l - 3a)3 ((1 - 9a)/s: + 6acl)Y^^ 
2 Ci,K.. ((1 - 9«)/«' + 2(k - 1 + 3a)c|,) V,' 



1 - 3a) (1 - 3a)3 

k((1 -9a)K + c|,(l + 3a)) 



a2 



((1 - 3a)3 ((1 - 9a)/t2 + QancD) 



2 A A 1/2 



cr 



(3.22) 
(3.23) 
(3.24) 

(3.25) 
(3.26) 



where Vs = ^iVj and V^ss = ^i^i^jj are gradients of the potential projected along the 
entropy direction in field space (and similarly in the following). The resulting equations 
of motion read, in Fourier space: 



z 
z 
w' 



0, 



,2,2 



+ a"/i; I Vs iva + c^^k v^ = Q 

w z 



(3.27) 

(3.28) 



where one can check that standard multifield DBI inflation [35] is recovered in the limit 
a = 0, K = 1. Similarly to the single-field case, one can also understand in a simple fashion 
the other extreme limit, in which one formally considers Mp = 0, by considering g as the 
cosmological metric. The proof, more subtle than in the single field-case, can be found in 
appendix (D). 

Let us now analyze the physical consequences of the above result. First, by taking the 
square roots in (3.19)-(3.20), we implicitly consider the regime in which the fluctuations 
are not ghosts, which reads 

a < (9 - Qc^/k)-^ . (3.29) 

In particular, a < 1/3 in the ghost-free regime, as we said in our discussion of the back- 
ground evolution in section (2.2). Interestingly, using its definition (2.27), the condition 



15 



(3.29) on a can be reformulated, through the first Friedmann equation (2.21), as an upper 
bound on the energy density pbrane, ^-e. there exists a critical background energy density 



above which the cosmological fluctuations become ghosts. This novel aspect due to the 
induced gravity can be simply understood in the formal limit Mp — )■ 0. Indeed, the bound 
(3.30) then boils down to V < 1//, which is a simple consequence of the form of the 
kinetic term in (3.17). 

Another point worth stressing is the very simple form of the coupling between adiabatic 
and entropy perturbations on super-horizon scales, i.e. when spatial gradients can be 
neglected: it is given by a single parameter ^ (3.26). Note also that this property is 
not particular to DBI Galileon inflation but is shared by a very large class of multifield 
inflationary models [116]. One novel feature however, due to the common presence of 
both the Einstein-Hilbert and induced gravity action, is the gradient coupling between 
adiabatic and entropy perturbations through the term c^^ in (3.27)-(3.28). When it can 
not be neglected, it implies that one must diagonalize the speed of sound matrix and that 
the adiabatic and entropy perturbations can not be quantized independently under the 
horizon (nor the field fluctuations themselves). In the following, we consider a regime in 
which this coupling can be neglected, in which case and are truly the respective speed 
of propagation squared of adiabatic and entropy fluctuations. Because s = cx)/{Hcx>), 
which can have a priori any sign and any value, enters into their explicit expression (3.22) 
and (3.23), and have no definite signs in general. With a slow- varying approximation 
in mind, we neglect s in the following, in which case one can check that and are 
both positive definite in the ghost-free regime, i.e. there is no gradient instability. In this 
regime, the entropic speed of sound Cg is always less than the adiabatic one c^, the only 
particular cases when they become equal being: 

• standard multifield DBI inflation (M — )■ and hence a — )■ and k — )■ 1). One then 
have = Cg = cd [34]. 

• the non-relativistic limit of DBI Galileon {cv — 1 for any M). One then have 

C(j Cg 1. 

• the limit in which the induced gravity dominates (Mj, -C M'^/cv and hence k ~ c^). 
One then have = Cs — 1 for any cv- 

Note also that ~ c^/k, is always less than one, contrary to the adiabatic speed of 
sound. If we require the latter to be subluminal, one should impose the condition a < 
(9 — 2c|,//t)~^, which is more stringent than the condition (3.29) for avoiding ghosts. 




(3.30) 
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Finally, by comparing the exact second-order action (3.21) with the one obtained from the 
procedure consisting in only perturbing the scalar fields, and not the metric: 



S, 



(2), naive 



^ (1 - 3a(3 - 24)) - 
1 - 3a / AT 2 i^Qs 



c 



(1 - «(5 - 24)) 



Q 



+ ... 



(3.31) 



we see that the latter is too naive and leads to completely wrong results (for example for 
the speeds of propagation) as soon as k differs significantly from one. This is in agreement 
with the structure of the solutions to the constraints equations in the flat gauge in appendix 
(B). 



3.1.3. Gravitational waves 

Here, we derive the second-order action for the tensor modes E. 



iJ5 



I.e. 



the transverse 



d'^Eij = - and trace free - i?* = - part of hij/a"^, where hij is the spatial metric in 
(3.1). This is actually straightforward within the bimetric approach. We know the result 
corresponding to the standard Einstein-Hilbert action: 



{2),EH 



dt d^k 



E 



2 ^2 



(3.32) 



where 



Ei, 



d^k 

(27r)3 



ikx 



(3.33) 



is the standard decomposition of Eij into + and x polarization states. Now, as far as 
the tensor sector is concerned, the Einstein-Hilbert and induced gravity and action are 
equivalent, with the replacements Mp -H- M and dt ^ dt, where t is the cosmic time 
associated to the induced metric, such that dt = cvdt (2.19). Hence, one immediately 
deduces from this the total second-order action for the tensor modes 



8 




Ml + M\v) {Elf 



(3.34) 



(this result actually also holds for a general, non constant, warp factor). In terms of the 
canonically normalized flelds in conformal time 



k 1 



(3.35) 
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this reads 



p 




(3.36) 



where q = ay Mp + M"^ /cv and 



Ml + M^cp 



(3.37) 



is the speed of propagation squared of the gravitational waves, which is always subluminal. 

3.2. Quantization and power spectra 

The exact calculation of the second-order action above made it possible to identify the 
remarkably simple form of the linear equations of motion (3.27)-(3.28). In the remaining 
of this paper, we restrict ourselves to considering a quasi de-Sitter inflationary phase - 
a(r) ^ — - in a slow- varying regime in which the time evolution of every quantity is slow 
with respect to that of the scale factor, such that for example z" / z ~ w" jw ~ ({' jq — . 
As explained above, we neglect the gradient coupling between adiabatic and entropy per- 
turbations, and we assume that the effect of the coupling ^ can be neglected on sub-horizon 
scales, so that adiabatic and entropy perturbations can then be quantized independently. 
We also consider an effectively light entropy perturbation, such that ^ 1, and 

which therefore gets amplified at (its) sound-horizon crossing^. 

The solutions to Eqs. (3.27) and (3.28) corresponding to the usual vacuum on very 
small scales then read 



form which one deduces, using (3.19)-(3.20) and (3.22)-(3.23), the power spectra of the 
initial adiabatic and entropy fields: 




(3.38) 



(3.39) 




(3.40) 



(3.41) 



Note that multifield models of inflation with intermediate masses of order H can have a rich phe- 
nomenology [117, 118], especially at the non- linear level. 
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Strictly speaking, the power spectra in Eqs. (3.40)-(3.41) are evaluated respectively soon 
after adiabatic and entropic sound-horizon crossing. Here, we have assumed that the 
hierarchy between the two sound speeds is not too large, so that the parameters entering 
into Eqs. (3.40)-(3.40), such as H and a, do not vary too much between the two epochs. 
Note that this is not as restrictive as it may seem, as the number of efolds between them 
depends only logarithmically on the ratio Cujcs- From the above results, one deduces that 
the enhancement of the entropy fluctuations with respect to the adiabatic ones that was 
found in standard multifield DBI infiation [34] - for which Qg ~ Qa/cv ^ is even more 
important in multifield DBI Galileon infiation, namely Vq* > ^Vq^. Finally, rewriting 
the above result (3.40) in terms of the curvature perturbation 

C = --Q., (3.42) 
a 

one finds its power-spectrum around horizon crossing 



V2vr; {{l-5a)K + 2aclf^ 

In the relativistic regime of particular interest to us, we have seen in section (2.2) that a 
phase of quasi de-Sitter expansion is achievable under the conditions (2.25) and (2.26). In 
this regime, the power spectrum (3.44) takes the form 

H V (l-3a) 



n.c,«i = ^^ 7^ (3-44) 

87r^Co-e \Mp J (1 — ba) 

The determination of the tensor power spectrum runs analogously. In the slow-varying 
regime, the Bunch-Davies vacuum solution of Eq. (3.36) reads 



from which one deduces, together with Eq. (3.35): 

2H^ {Ml + M'^/cv 



,1/2 



{Ml + M^cvfl-" 



In the relativistic regime, in which the condition ^ Mlc^ should hold to achieve 
an infiationary phase, the effect of the induced gravity is therefore subdominant, and the 
power-spectrum of the gravitational waves is not modified at leading order compared to 
its result in Einstein gravity. 
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3.3. Large scale evolution 

In the single field case, relating the quantum fluctuations generated during inflation to 
the primordial seeds of the large-scale structure deep in the radiation era relies on the 
constancy of the curvature perturbation ( on super horizon scales. As is well known, and 
shown for the flrst time in [119], this property does not hold in general in a multifleld 
situation due to the large scale feeding of the curvature perturbation by the entropic 
fluctuations. This can easily be seen here by noticing the existence of an exact flrst-order 
integral of the adiabatic equation (3.27) when spatial gradients can be neglected^: 

^-v, ^ C ^ --Vs . (3.47) 
z / z z 

As for the entropic perturbation, one can use this result, together with Eq. (3.28), to show 
that it evolves independently of the curvature perturbation on large scales: 



w" 



< + (^-— + +ejvs^O. (3.48) 

Quite generally, upon neglecting the decaying modes, one can therefore express the late- 
time curvature perturbation as 

C — '^aQa* ~l~ ^sQs* (3.49) 

where 

UA ^^"'^^l a ((l-5«).^ + 2a.4)3/4 )^ ^^"^"^ 

and the complicated expression of As has been chosen such that 

Pc = n. , (3.51) 

i.e. measures the contribution of the entropy perturbations to the flnal curvature 
perturbation. In the following, we concentrate on the curvature perturbation and leaves 
aside the interesting possibility of a remaining primordial isocurvature perturbation. 



4. Non-gaussianities 

This section is devoted to the analysis of the primordial bispectrum generated by mul- 
tifleld DBl Galileon inflation in the relativistic regime. Schematically, primordial non- 
Gaussianities can be generated at three different epochs of the inflationary evolution: 



^ This result can also be obtained from the large-scale limit of the energy constraint (B.8) together with 
the momentum constraint (B.7). 
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• before the relevant modes cross their respective sound horizon. The fluctuations are 
then in their vacuum and their rapid oscillations usually average out. However, in 
cases in which the background evolution displays sharp or periodic features, construc- 
tive interferences between different modes can be generated, and hence significant 
non-Gaussianities [120-124]. 

• around horizon crossing. This contribution is unavoidable and is typically large 
in models in which derivative operators are important, like in scenarios with non- 
standard kinetic terms (see e.g the review [125]). 

• after horizon crossing. The curvature perturbation being non-perturbatively con- 
stant on large scales in single-field models [113, 114, 126-130], this mechanism is 
only possible in multifield scenarios. The resulting shape of the non-Gaussianities 
depends on the masses of the entropy fluctuations while its amplitude is very model- 
dependent (see e.g [131] for a recent illustration). 

Following our analysis of the linear fluctuations, we consider a smooth quasi de-Sitter 
background evolution and hence do not consider the first type of non-Gaussian contribu- 
tions. As for the ones after horizon crossing, they would be of the well known local type 
as we consider en effectively light entropy perturbation. Their amplitude being model- 
dependent, we do not consider them in the following and concentrate on the second type 
of contributions (Ref. [45] gives an example, in the DBI context, in which local- type 
non-Gaussianities and ones generated around horizon-crossing can both be significant). 



4.1. Third-order action 

To determine the bispectrum generated by the quantum interferences amongst different 
modes of perturbations, we should calculate the third-order action in these perturbations. 
Fortunately, in the almost de-Sitter and relativistic regime of interest (see section 2.2), we 
do not need its exact form. Indeed, using the results of appendix (B), one can check that, 
in the flat gauge, the metric perturbations are negligible compared to the field fluctuations 
themselves. Therefore, it is sufficient for our purpose to perturb the scalar fields only and 
at leading order in the slow-varying approximation. Hence, one can consider separately 
the brane action (2.5) and the induced gravity action. The corresponding result for the 
former, determined in [35], is given by 



+<^Q.Ql + <^q}^^ - 2^^dQ^dQsQs ] . (4.1) 



For the induced gravity action, it proves very useful, like for the second-order action, to 
start from its ADM form Eq. (3.6). Here, we simply give the result, leaving details of the 
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calculation in appendix (C): 



S'^^l = [ dtd'xa'fH'M'a f-^(5 - 2cl)Ql + ^-^^qJ 



a2 



6 ^2d''Q. 3 {dQ^fd^Q^ 



X) 



For completeness, we have kept the terms subleading in c|, that appear in the calculation 
(in the first line). However, we should bear in mind that terms not calculated here, coming 
from the mixing with gravity or beyond leading order in the slow-varying approximation, 
can be of the same order of magnitude or greater as those terms (see for instance [132, 133]). 

Note that, at leading order in a slow varying approximation, only shift-symmetric 
operators can appear, and only with an even number of entropic fields Qs- Note also that 
the seven mixed adiabatic/entropic operators are the multifield generalisations of the four 
purely adiabatic ones. Eventually, one can check that the relation between the adiabatic 
speed of sound and the coefficient of the operator Qa{dQa)^ that was noted in the 
context of the effective field of theory of infiation [134] is verified. However, because of 
the presence of higher-order vertices, the direct link between and the amplitude of the 
equilateral-type non-Gaussianities is spoiled, as has already been noted in [67]. To see 
this, let us use the (leading-order) equations of motion 



2 



Q^ + 3HQ^-ci—f^ = (4.3) 



Q^ + 3HQ,-cl^ = (4.4) 



-2 _ 

,10 



to simplify the third-order action . Using integrations by part, one finds that 



dtd'xa-'Qi^P = j dtd'xa' (^^^-J (^•^) 
y dtd3xa3-__ = y dtd3xa3^-g3 + _g^_j . (4.6) 

Hence, one can trade the two higher dimension operators Q^d^Qa and {dQaYd^Qa in the 
second line of (4.2) for the standard fc-infiationary operators and QaidQuY-, as was 



The fact that this procedure enables one to calculate accurately correlation functions (if boundary terms 
are carefully taken into account) was discussed in Refs. [135-137], more recently in Refs. [133, 138, 139] 
and applied to the most general scalar- tensor theory with second order equations of motion in [140]. 
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noted in Ref. [70]. Similarly, for the mixed adiabatic/entropic operators, one can show 
that 



„3 A A '^^Qs 



dt d-^a; cl^ Q^Q 



2 



\ c2 



and 



/ 



1^13 3 1 dQadQs d'^Qs 
at a x a 



H a? 



dt d^x a^- 



2H^' a? H-^" -2 



H 



2H a2 a2 



(4.7) 



where we have introduced the ratio between the entropic and the adiabatic speed of sound 



Co- 



'1 -3a(3-24) 



1 - a(5 - 2c 



(4.9) 



V) 



Using the four equalities (4.5)-(4.8), one can then write down the effective total third-order 
action as 



rieff 
•^(3) 



, „4/. (^^l)!_24(l-2a)^^4 



2cj) 



— -ac 



■D 



1 



1 Q; 



A2 



1 



a? Ha? 



(4.10) 



with 



= 1 - 3a (5 - 24 - 4A^ + A^) 
^Q.iSQ.)^ = 1 - a (9 - 24 - 3A2) . 



(4.11) 
(4.12) 



As we will see, the correlation functions induced by each of the seven operators in (4.10) 
are linearly independent. Hence, no integration by part can simplify this result further, 
i.e. none of the operators is redundant [140]. In particular, contrary to the purely adia- 
batic perturbations, one can not get rid of all the higher-order operators in the multifield 
scenario, hence the last line in Eq. (4.10). 

At lowest order and at tree level, one deduces the bispectrum from the third-order 
action as [141] 

(C(t, fci)C(t, fc2)C(t, k,)) = -I f dt'(0|C'(t, fci)C^(t, k,)C\t, k,)Hj{t')\0) + c.c (4.13) 

J —00 
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where is in the interaction picture, Hj denotes the interaction Hamiltonian, such that 
S^^^ = - J dt'Hi{t'), and |0) is the free theory vacuum. In practice, one evaluates this 
integral using the conformal time r ^ —l/[aH), in which case one can extrapolate the 
upper bound to r = as the main contribution to the integrand comes from the period 
around horizon crossing. One also uses the appropriate rotated contour in the complex 
plane (r — )■ — (oo — ie)). Following this procedure, and using Eq. (3.49) and the structure 
of S^^y the bispectrum is calculated as 

(C(fci)C(fc2)C(fc3)) = {A^)''{QAki)QAk2)QAk3)) 

+ AMs?{{QAki)Qs{k2)Qs{k3)) + perm.) , (4.14) 

where the 'perm.' indicate two other terms with permutations of indices 1, 2 and 3, and 
whose result we conventionally write as [125] 



i=l 



3) 



where is the primordial power spectrum given in Eq. (3.51) and the dimensionless 
shape function 

has been separated out between its purely adiabatic and entropic-induced components, 
corresponding respectively to the first and the second line in (4.14). A peculiarity of 
standard multifield DBI inflation (corresponding to = 0) is that 5''^°"'^'^^ = 5''^'^''**^. 
However, as we will see, this is not the case in multifield DBI Galileon inflation. Hence, we 
analyse 5*'-'^°''^^ and S^'^'^^^ separately in the following, considering first the purely adiabatic 
signal and then the entropic-induced one. 

4.2. The adiabatic bispectrum 

The adiabatic shape 5*'-'^'^'^^ is given by 

5'('^'^'^ = SAqsSqs + ^^Q^(aQ^)25'Q^(5Q^-)2 (4.17) 



where Aqs and ^Q^(aQ^)2 are given in Eqs. (4.11) and (4.12) respectively and 

5q3 = Si{h,k2,h) (4.18) 
^Q.idQ.r = '^2(A;i, k2, h) + 2 perm. (4.19) 
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with (this will be useful in our discussion of the entropic-induced shape) 

Si{h,h,h) = '^ (4.20) 
S2{h, h, h) = ■ k2 {2hk2 - k,K + , (4.21) 

and K = ki + k2 + k^. The adiabatic shape is thus a linear combination of the well known 
shapes of /c-inflationary type Sqs and Sq^^^qq^^, which are displayed in Fig. 1(a). One can 
see that there are very similar, in particular both peaking on equilateral triangles ki = 
^2 = k^. To assess the experimental ability to distinguish between different bispectrum 
momentum dependence, we use the scalar product between shapes introduced in [142] (see 
[143] for earlier similar definitions). In this sense, the correlation C[Sqs, Sq^(^qq^'^2] equals 
0.97, which explains why the two shapes are often approximated in CMB data analyses 
by a common (factorizable) ansatz, called equilateral [144]: 

= - (^Yk + 2 Perm.^ + (f ^ ^ perm.^ - 2 . (4.22) 

Its shape is displayed in Fig. 1(b) (left), where one can see that it captures their main 
features. Indeed, C[Sq3 , Sec\ = 0.94 and C[Sq^(^qq^-^2, 5'eq] = 0.99. However, as pointed out 
in [132], Sq3 and Sq^^qq^^ are different, especially for squashed triangles ^2 + ^3 = ki/2 
and more generally for flattened triangles ^2 + ^3 = ^1. One can therefore highlight their 
differences by considering an appropriate linear combination of them with respect to which 
it is almost orthogonal. This was carried out in Ref. [145], in which the authors designed 
a corresponding factorizable template, called orthogonal: 

5orth = -3 + 2 perm.^ + ^ + ^ P""™') " ^ ' ^^'^^^ 

that is represented in Fig. 1(b) (right) (we actually plot its absolute value so that 
its difference with the equilateral ansatz for flattened triangles is more visible). This 
orthogonal shape is now routinely considered in CMB data analyses, which put con- 
straints on its amphtude f^l^ (defined such that S in Eq. (4.15) equals jQfNL^Sorth)- 
f^l^ = -79.4 ± 133.3 (95% C.L) (Fergusson et al. [3]) and /^''f = -202 ± 208 (95% C.L) 
(WMAP7 [4]) . 

However, a concrete early universe model that generates such a non-gaussian signal 
with an observably large amplitude was still lacking^^. We now demonstrate that DBI 



Shapes referred to as orthogonal in Ref. [146] reach their maximum in the squashed configuration, like 
the orthogonal shape. However, their sign is the same for all triangle configurations and their correlation 
is large not with the orthogonal template but with the equilateral one. 



25 




(b)Shape of 5cq in Eq. (4.22) (left) and of the absolute value of S'orth in Eq. (4.23) (right). 



FIG. 1: Different shape Sx as a function of We set tliem to zero outside tlie 

region 1 — k2/ki < k^/ki < ^2/^1 and normalize them to one for equilateral triangles ^ = 

Galileon inflation provides such a model. The reason is simple: as a increases, Aq3 and 
^Qa(9Q<T)2 deviate from 1 but, whereas the latter remains positive and tends to 0, the latter 
becomes negative (see Fig. 2). Hence, the adiabatic shape (4.17), positively correlated 
with the equilateral template when a = 0, becomes anti-correlated with it as a —t- 1/9 
(see Fig. 3). As a consequence, in a transient region - centered around a = 0.097 - 
C[S^"°''^\ Seq] becomes small and the similarities between the two equilateral-type shapes 



26 




0.00 0.02 0.04 0.06 0.08 0.10 



FIG. 2: Aq3 in Eq. (4.11) (dashed red) and ^q^(^qq^)2 in Eq. (4.12) (plain blue) as we vary a 
with c|, <C 1. 



are efficiently subtracted, resulting in a large correlation with the orthogonal template of 
maximum value 91 %. One should note that this value is not surprising, as it is the value 
already given in Ref. [145] for the correlation between the true orthogonal shape - arising 
as linear combinations of Sqs and Sq^i^qq^^ - and its factorizable template Eq. (4.23). 
The corresponding shape is represented in Fig. 4 (left). 



As far the amplitude of this non-Gaussian signal is concerned, we find that^^ jorth [aacr) 



' NL 



0.016 



(at a = 0.097). Outside the transient region - for definiteness let us say for 
a < 0.095 and a > 0.100, the adiabatic shape can be considered as truly equilateral, and 
its amplitude simply measured by evaluating ^S" on equilateral triangles. One then finds 



neq (aaa) 
J'Nh 



1 



(21 - 404a + 2233a2 _ 30660=^ 



324c2 (1 + T2)2 



[1 -5a)^(l -9a) 



(4.24) 



whose dependence on a is represented in Fig. 5 (left). Consistently with the correlation 
of the adiabatic shape with the equilateral ansatz (top of Fig. 3), one finds that, as a 



increases, /] 



eg (crcrcr) 
NL 



goes from negative to positive values. It also grows unboundedly as 



a — >■ 1/9, corresponding to the fact that the kinetic term of the adiabatic fluctuation 
becomes close to zero and that the theory becomes strongly coupled in this regime. To be 
more quantitative, the figure 6 (left) gives the region allowed by the current observational 
bounds on non-Gaussianities in the two-dimensional parameter space [cd, a) for T^g = 0. 
We used the constraints f^l = 143.5 ± 151.2 and f^''^^ = -79.4 ± 133.3 (95%C.L) from 
Fergusson and Shellard [3] and took into account the condition that the shape is sufficiently 
correlated (or anti-correlated) at 80% with the equilateral or orthogonal templates. The 



This value is calculated as f'^l'^ — ^ c[S^'^tf°s ''L] ' ^ relevant to compare to the observational 
constraints [3, 4]. Alternatively, one can estimate it by simply evaluating on equilateral triangles. 

0.011 

'cl,(l+T2J2- 



One then finds /^''f ^"""^ 
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FIG. 3: Correlation of the adiabatic shape S^"''^'^') with the equilateral ansatz (top) and the 
orthogonal ansatz (bottom) as a function of a in the relativistic limit c|, <C 1. 



unbounded growth of f^\}'^'^'^^ with a of course explains the sharpness of the constraints 
when a > 0.10. Note also that smaller values of cv are allowed for a ~ 0.097 compara- 
tively to other values of a because of the small numerical factor in f^l^ = — tjutft-^^ 
which itself is explained by the partial cancellations between Sq^ and Sq^^^qq^^ required 
to generate orthogonal non-Gaussianities. 

We have concentrated so far on the equilateral and orthogonal shapes of non- 
Gaussianities. However, one should point out that shapes yet qualitatively different from 
them do arise in our model. In particular, as is clear from fig. 3, neither the equilateral 
nor the orthogonal templates faithfully represent the adiabatic shape at a ^ 0.0975 (the 
corresponding correlations being respectively —0.55 and 0.74). Indeed, the corresponding 
shape, represented in Fig. 7, is very close to what the authors of Ref. [145] have called the 
flat shape. To be complete, one should actually note that shapes very similar were shown 
to arise in different contexts: the "surfing shape" of Ref. [70] or the "flat" shapes gener- 
ated by curvature-related terms in [147]. Eventually, all of those shapes also have a large 
overlap with the enfolded template [148], aimed at modeling the effect of a non-Bunch 
Davies vacuum [132, 149], and which is actually a linear combination of the equilateral 
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FIG. 4: Absolute value of the adiabatic shape S^'^"'^^ for a = 0.097 (left) and of the entropic 
shape S^'^^^^ for a = 0.080. We use the same conventions as in Fig. 1. 




/n: 




FIG. 5: /n'l '""'^ in Eq. (4.24) as a function of a for cv = 0.1 and T^^ = (left), f^l^"''^ in Eq. 
(4.35) as a function of a for cx> = 0.1 and T^g = 1 (right). 

and orthogonal template [145]. For instance, our shape in Fig. 7 is 93 % correlated with 
it. 



4.3. The entropic-induced bispectrum 

The entropic-induced shape 5'^'^'^*'' in Eq. (4.16) reads 



)2 

,)2 + - '^Ol)SdQ^dQsQs 



a 



3 A Sq2q2q^ - S(^9Q^)2g2Q^ 

(4.25) 
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FIG. 6: The region in the parameter space (cD,a) allowed by the observational constraints on 
f^l and f^^l^ for T,, = (left) and T,, = 10 (right). 



with 
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X 



3a 



S. 



s 



dQa-dQsQs 



1 - 3«(3 - 2c|))(l - «(5 - 24)) 



— S'i(A;i, AA;2, \h) + 2perm. 
82(^^2, Xk^, ki) + 2 perm. 

\s2{ki, Xk2, Afcs) + Sperm. 
A 

kik2k^ 



Kf 



{Ki + 3ki) + 2 perm. 



(4.26) 

(4.27) 
(4.28) 

(4.29) 
(4.30) 



S, 



ki fca ■ fcg 

k2k, Kf 



[Kf + kiKl - {kl - \^k2k^)Ki + 3X^kik2k3) + 2 perm. (4.31) 



with Ki = ki + A(/c2 + ^3). The fact that x verifies < x < 1 with x = 1 when a = 
imphes that the induced gravity tends to lower the intrinsic ampUtude of quantum non- 
Gaussianities induced by the entropic perturbations compared to the ones generated by the 
adiabatic perturbations. However, depending on the magnitude of the transfer function 
Tg-s; and hence on the large-scale evolution, the total bispectrum can well originate mostly 
from entropic origin. 
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FIG. 7: The adiabatic shape S^"'^'^^ for a = 0.0975. We use the same conventions as in Fig. 1. 

Similarly to the purely adiabatic sector, in the relativistic regime c|, ^ 1, the entropic- 
induced shape 5'^'^'^'*^ represents a-one parameter family of shapes depending on a. How- 
ever, it is also interesting to study each of the shapes (4.27)-(4.31) separately^^. One should 
note then that one can not factor out A from their fcj's dependence. In other words, these 
shapes constitute themselves five families of A-dependent shapes, which one can check to 
be linearly independent. In this respect, A = 1, corresponding to standard multifield DBI 
inflation in which the adiabatic and the entropic speed of sound coincide, is an enhanced 
symmetry point at which the entropic-induced shapes can be expressed in terms of the 
two purely adiabatic ones (4.18) and (4.19): 



A 



-S. 



6 ^ 

Is 

2'^dQ^dQsQs 



= Sr 



(4.32) 

(4.33) 
(4.34) 



However, for general values of A, some of the entropic-induced shapes differ significantly 
from the equilateral template. The most dramatic changes occur for Sq^q2 and Sq2q2q^, 
which are represented in Fig. 8 for A = 0.1. One can see that these shapes reach their 



There is some arbitrariness here because, as we have pointed out, some vertices, and hence some shapes, 
may well have been replaced by others by using the linear equations of motion and integrating by parts. 
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FIG. 8: Shape of Sq^q2 in Eq. (4.27) (left) and Sq2q2q^ in Eq. (4.30) (right) for A = 0.1. We 
use the same conventions as in Fig. 1. 

maximum near the squeezed limit fcs <^ A;2 — ki (all the more so as A decreases), which is 
unusual for non-Gaussianities generated around horizon crossing^^. Interestingly, note that 
the authors of Ref. [133] have recently shown that similar shapes (referred to as S'^ and Ss 
in their paper) arise at next-to-leading order in slow-roll in fc-inflation. Now coming back 
to the global entropic-induced shape S^'^^^\ one can see from its expression (4.25) that the 
contributions from the shapes Sq^q2 and Sq2q2q^ are diminished by factors of A^ compared 
to others. Hence, we do not expect the total entropic-induced shape S*^""***-* to exhibit 
features similar to the ones in Fig. 8. Indeed, representing it graphically shows that it is 
qualitatively similar to the adiabatic one: it interpolates between negative (when a = 0) 
and positive equilateral non-Gaussianities passing through orthogonal ones, although for 
smaller values of a centered around 0.08 (the corresponding shape is represented on the 
right of Fig. 4). This is confirmed by the correlations of S**-""*^-* with the equilateral and 
orthogonal templates, shown as a function of a in Fig. 9. Let us note also that S^'^^^^ is 
of "surfing" type, similar to the adiabatic shape in Fig. 7, for values of a slightly larger 
than the transition value 0.08. 

Regarding the amplitude of the entropic-induced non-Gaussianities, we find that 
j.mi,h{ass) _ _ 0-0043 for a = 0.08 while it can be simply characterized by evalu- 



One should note however that, as expected, they vanish as k^/ki in the squeezed limit, like all the other 
shapes generated in this model. 
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FIG. 9: Correlation of the entropic-induced shape S^"^^^ with the equilateral ansatz (top) and 
the orthogonal ansatz (bottom) as a function of a in the relativistic limit c|, <C 1. 

ating yS" on equilateral triangles for a < 0.078 and a > 0.082, leading to 

re,iass) ^ 5 v/(l-9a)(l-5«) 

/nl ^ y2j2 12(1 + 2X)^cl (1 - 3a) (1 - 9a) (1 - 5a)2 

(33 + 30A - a(762 + 648A) + 3a2(1631 + 1282A) - a^(6532 + 5276A)) , (4.35) 

where A ~ a/(1 — 9a)/(l — 5a) in the relativistic regime. The dependence of f^^'^'^^^ 
on a, similar to the one of f^j^'^'^"^ in Eq. (4.24), is displayed in Fig. 5 (right). The 
observational constraints on the parameter space (cd, a) for To-s = 10 - in which case 
the total bispectrum is dominated by its entropic-induced component - can be seen in 
Fig. 6 (right). Values of cv smaller than in the case of purely adiabatic non-Gaussianities 
are allowed because of the suppression factor ^ but the pattern of constraints is 

qualitatively similar. 

In all our above discussion on the non-Gaussianities generated by multifield DBI 
Galileon inflation, we considered the relativistic regime c|i -C 1 in which large non- 
Gaussianities are generated. As we have explained, for bigger values if c|,, the gravitational 
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back-reaction can not be neglected and our analysis is not complete. However, to get a 
taste of the shape of non-Gaussianities for intermediate values of c|,, it is interesting to 
take the non-relativistic limit cx) — ?■ 1 (and hence A ~ 1) in our expression (4.16) for the 
shape. The overall amplitude of non-Gaussianities become of course small in this limit (as 
1 — c|,) but interestingly, the corresponding shape simplifies to 

^^T^ = -2(l-3a)(l+V,^J ('^«^ + - • ^^-^'^ 

In particular, the entropic-induced shape and the adiabatic shape become equal in the 
non-relativistic limit, so that the total shape becomes independent of the transfer function 
To-s, as in standard multifield DBI inflation. Note also that, as | — 2a > 0, there is no 
cancellation between 5*^3 and Sq^(^qq^y in Eq. (4.36), so that this shape is of equilateral 
type for any values of the induced gravity strength a. Of course, these results should be 
taken with a grain of salt as other terms that we have neglected compete in magnitude 
with the result (4.36) above. Nonetheless, they give an indication that the rich structure 
of shapes that we have studied in this paper does not persist in the non-relativistic regime. 



5. Conclusions 

In this paper we have studied the cosmological implications in the early universe of the 
multifield relativistic Galileon model in which an induced gravity term is added to the 
DBI action governing the motion of a D3-brane in a higher-dimensional spacetime. By 
employing a bimetric perspective in which the induced gravity and Einstein-Hilbert action 
are treated on equal footing, we have been able to analyze this model in a physically 
transparent and computationally very efficient way. After deriving the gravitational and 
field equations of motion in full generality, we have considered the background homogenous 
evolution. We have shown that the induced gravity tends to violate the null energy 
condition and to hamper infiation in the relativistic regime in which the brane almost 
saturates its speed limit. Although the induced gravity should thus be subdominant for 
our model to sustain a phase of quasi de-Sitter expansion in the relativistic regime, we 
have shown that its effects need not be negligible. 

Our study of the cosmological perturbations generated during infiation even revealed 
that the induced gravity can leave very non-trivial imprints on the cosmological fiuc- 
tuations. We computed the exact second-order action in the fiuctuations, showing the 
existence of a critical background energy density above which cosmological fluctuations 
become ghosts. We have also shown how the fiuctuations react both to the cosmological 
and to the induced metric as a consequence of the bimetric structure of our model. Re- 
stricting our attention to the case of two fields, we were able to cast the exact coupled 
linear equations of motion for the adiabatic and entropic perturbations in a very elegant 
and intuitive form. An interesting feature of our model lies in the fact that neither these 
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perturbations not the field fluctuations themselves are decoupled on sub- horizon scales. 
Leaving the study of its consequences to future works, we neglected this coupling and 
derived the properties of the cosmological fluctuations at leading order in a slow-varying 
approximation. Wc then showed that the cntropic speed of sound differ from (and is less 
than) the adiabatic one. As a consequence, the two types of perturbations are amplified 
at their respective sound-horizon crossing and the amplitude of entropic perturbations are 
even more enhanced, compared to standard multifield DBI infiation, than the adiabatic 
ones. 

We then considered the non-Gaussianitics of the curvature perturbation generated in 
the rclativistic regime. For that purpose, we computed the dominant contributions to the 
third-order action and we identified its seven independent cubic operators: two purely 
adiabatic and five mixed adiabatic/entropic ones. As a consequence, the shape of the 
primordial bispectrum induced by the entropic fiuctuations is different from the purely 
adiabatic shape, thereby breaking the degeneracy of standard multifield DBI infiation. 
Each of these components of the bispectrum depends non-trivially on the strength of the 
induced gravity as measured by a parameter o; < 1/9. Although quantitatively different, 
they share common features: their shape is of equilateral-type for small and large (the 
largest allowed) values of a, f^j^ being then respectively negative and positive. In transient 
regions, centered respectively around a — 0.097 and a — 0.08, the two shapes become of 
orthogonal type. We stress that it is the first concrete model in which this shape of 
non-Gaussianities, which is routinely considered in GMB data analyses, is generated with 
an observably large amplitude. More generally, we find quite remarkable the fact that 
qualitative aspects of the bispectrum, such as its shape and sign, depend very sensitively on 
the precise value of the induced gravity strength. This exemplifies the usefulness of higher- 
order statistics for pinning down the mechanism that seeded the large-scale structure of 
the Universe. 

Eventually, the fact that the induced gravity tends to violate the null energy condition 
is potentially interesting in the context of bouncing cosmologies, in which this condition 

should be broken in order to connect a contracting to an expanding universe. We leave 
for future studies the interesting question of whether a stable bounce solution can be 
constructed in DBI Galileon models. 
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Appendix A: Details on the explicit form of the action 

In this appendix, we give the exphcit form of the induced gravity action S'ind = 
^- J d'^Xy/—'jR['j] in terms of the scalar fields and the geometrical quantities associated 
to the cosmological metric, leading to a multifield relativistic extension of the quartic 
Galileon Lagrangian in curved spacetime. 

Remembering that 'j^u = h^^/'^g^y where g^y = Qnu+f Gijd^cj)^ dy(t)'\ one readily obtains 

Sind = J ^ x — ^/-gR[g\ - J d x^-g— J^^g (^-1) 

where Xl'^ = —\g^''d^(i)^d^(i)'^ and = M'^/\/h. The real difficulty is to obtain an 
explicit expression of the inverse g'^^ of g^i, as well of its Ricci scalar, the two metrics g^^ 
and g^j^u being related by a disformal transformation. Using the results in Ref. [150], it 
turns out one can obtain the following non-trivial expression for g^'^: 

-gt^u = Agf^'^ - AjV^0^V'^0'^ (A.2) 

where 

A ^ i (l - 2/X/ + AfX^lx'} - S/^xf Xixjl) (A.3) 

and 

Au^^[[l- 2fXi + 4/2xlf Xgj) Gij + 2/ (1 - 2/Xf ) Xjj + AfXjKXf^ . 

(A.4) 

Here and in the following, field space indices I,J... are lowered (respectively raised) with 
the field space metric Gjj (respectively its inverse). V is given explicitly in Eq. (2.9) in 
terms of the kinetic terms X^"^ defined in Eq. (2.8) and denotes the covariant derivative 
associated to the cosmological metric g^^^. Note also that all this complexity disappears 
in the single field case, in which A ^ 1, Au — )• //(I — 2fX) and D — )• 1 — 2fX. 

To calculate R[g], one possibility is to use the formulas relating different covariant 
derivatives and their related geometrical quantities, as given for example in Ref. [151]. 
Here, we rather follow a more down-to-Earth approach. We start from the general expres- 
sion 

Rxi^uM = \ i9f.u,KX + gxK,^,u - gxu,^f. - 9^.K,ux) + 9r,a{A9\lxA9]1u - AaZxAaMn) (A-S) 

that we evaluate in a locally inertial frame (in which d^g^p = and hence r[5f]^^ = 0) 
before covariantising it. In such a frame, one finds that 

= fiAGjj + 2^,^GjMX^^)V>^(a^a>-^ + fi^Va^^V^^^) (A.6) 
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where a hat denotes a quantity defined with respect to the metric Gjj = fGjj, so that 
1 1 

^AB ^ '^{^Af,B + ^bIa - G f,lGAB) + -jG {GcA,B + GcB,A- Gab,c) ■ (A.7) 

In this frame, one also has 

9ij,v,k\ ~I~ 9\k,ij,v 9\v,Kii 9lJ.K,v\ 9fiiy,K\ ~(" 9\k,^u 9Xu,kii 9hk,v\ 

+ 2RujK<p^jfx<l>',<f>-!u (A.8) 

and 

~9vAm:xmi. - mi.miAQ]) = -2(0^,. + r^.<[x'^fj (^Si. + ^cD^%]<f>^.){HAB - jgab) 

with 

Hab = JGab + '^A^ f'^{XAB — '^fXAcX^) + ^Af"^ Xi^AC'^^'^' {Xb)E — '^fXff^XEH) 

+ ^fA^'^A'^^XAcXBFiXEG - VXehX^) . (A.9) 

Covariantizing the result, i.e., with the substitution 0^^^^ — )■ V ^jSJ^cf)^ , one deduces the 
compact expression 

where 

= + fi^V^0^V.0^ . (A.ll) 

From this, one easily finds the explicit expressions for the induced Ricci tensor and Ricci 
scalar (we use the convention that for any rank-2 tensor 1^,^, \Y] = g^'^Y^i,): 

R^p[g] = AR^fsig] - AijV''(l>'V^<j>'Rxau0[9] 

+ ahab ([n^nf^ - (n^ • u^u) 

- 2{AX'^' + 2X'^''X'''AMN)RLijKVp(l>''Va(t>' (A.12) 

and 

R[g\ = A^R[g\ - 2AAij{d(l>' ■ R[g\ ■ dct)') + AijAMLV^4>'^''4>''^^(t>^^''(i>''Rx^.M 

+ a'Hab ([n^][n^]-[n-^-n^]) 

- 2AAcD Hab {[U^'W ■ • 90^) - (^c/.^ • • • a</.^)) 

+ {AijAkl - AjlAkj)Hbe {d(p' ■ • dct>'){d^'' ■ • 90^) 

+ AiAX"^' + 2X'^''X'''Amn){AX''' + 2X'''X''''Acd)Rlijk • (A.13) 
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Eventually, we have thus obtained the explicit expression 

Sind = / d^x {C-a + U + Cc + Cd + Ce + Cf + Cg + U) (A.14) 



where 



Ca = A^VVR[g] (A.15) 

£b = -2AVVAij{d(j)^ ■ R[g] ■ dcj)'^) (A.16) 

= a'VvHab ([n^][n^] - [n^ ■ n^]) (a.it) 

Cd = -2aVvAcd Hab ([n^](50^ ■ ■ -90^) - (^^^ ■ ■ ■ 90^)) (a.is) 

= Vv{AuAkl - AilAkj)Hbe (90^ ■ ■ 90^) (90^ • ■ 90^) (A. 19) 

Cf = v^^,j^mlVV'V^0^V^0^^V>^/?a^,«[(?] (A.20) 
Cg = 4yp(AX^^ + 2X^*'X-^^^M7v)(AX^^ + 2X'^X''''Acd)Rlijk (A.21) 



C, = {AX'' + ^AklX'^'x'^') (A.22) 

and Ce.Cf and Cg represent purely multifield effects. One should also note that pertur- 
bative calculations hint at a fully non-linear equality between Aab and Hab although we 
were not able to prove it in the general case. 



Appendix B: Solutions to the constraint equations 

In this section, we give details about the constraints equations (3.10) and (3.11), their 
solutions in different gauges, and the calculation of the second-order action. 

• Single field case and uniform infiaton gauge 

In the single field case, we write the scalar field Lagrangian as P(X, 0) where X = 
— \g'^^d^(j)di,(j) and work in the uniform infiaton gauge in which hij = a{tye'^^6ij and the 
field takes its unperturbed value. The linearised momentum constraint (3.11) then gives 
Vi = and 

5iVuni = |:«:, (B.l) 

where k was defined in Eq. (2.23). The linearised energy constraint (3.10) then yields 

/3uni = — ^K, + ax (B.2) 
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where 



Ml 



Ml 



M2 



c 



2 ^2 + "^l^ '^^ + A#2^3 



2 /7'2 /vf 2 „: 



(3-4)«: 



2r2 



(B.3) 



and l/c|_inf — 1 = 2XP^xx/ P,x- Expanding the action to quadratic order in C, by substi- 
tuting^^ the expression (B.l) for 5N (as usual, the exphcit expression of /3 is not needed 
for the second-order action), one obtains 



(2) 



()2 k'^P 



X 



2Hmi cti„f 



3(1 - + 



3M2 
Mlcv 



1 + 



-d + e) 1 + 



K — 



1 + 



1 - 
Mlcv 



(B.4) 



where we have used the two background equations 



and 



Ml + 



Ml + 
M^ 



M^ 



3H^ + P- Px0^ = 



a 

2- + H' 
a 



2H 



Cv Hcv 



-P. 



(B.5) 



(B.6) 



in particular to show that the terms in (^^ identically vanish. In our case of interest for 
which P = -1//(V1 -2/X -l)-V, this yields the resuh Eqs. (3.13)-(3.15). 

• Multifield case and fiat gauge 



In the multifield case and in the fiat gauge in which hij 
momentum constraint (3.11) gives Vi = and 



a{t)'^6ij , the linearised 



2HMl 1 + 



M^ 



clMl 



6Ni 



flat 



3a) -^rQ' + ^hQ' 
Cv cvH 



(B.7) 



while the linearised energy constraint (3.10) yields 



2H 



Mi + 



+ 



+ 5Ni 



flat 



a 



— 3" + QH'Ml I 1 + 
^v 



+ (1 - 9a) ^ + VjQ' 



2c%Ml 



0. 



(B.8) 



In the ADM formalism, it is sufficient to use the perturbed lapse and shift up to first order as their 
second-order parts cancel out in the action [107, 132]. 
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Note that, using the background equations of motion, one can check that the two above 
solutions agree, in the single field case, with the results obtained respectively from Eqs. 
(B.l) and (B.2) by making the explicit gauge transformation between the uniform infiaton 
gauge and the fiat gauge: 



6N, 



flat 



1 d 
Hdt 



H 



-Q 



11 



(B.9) 



and 



'flat 



acj) 



— a 



Ml 



X 



2H^Ml cl 



+ 3(l-/€) 



3M2 



1 - 



(3 - cI)k 



2r2 



d 




dt 


L0 J 



Eventually, let us give the form of the "anisotropy constraint", i.e. the traceless part of 
the linearised gravitationnal i — j equation of motion: 



Ml + —\[ 5Nfi,t + /3flat 



Cv 



ff/3flat ) + 2if/3flat I Ml + (1 - s/2) 

Cv 



1 



s \ fM^a. ^ 
^2 / e^/V H : c^jQ = . 



V 



Cv 



(B.IO) 



Given that the expressions in the fiat gauge of two gauge-invariant Bardeen potentials are 
given respectively by $ = 5A^flat + /3flat and \1/ = — -ff/3flat, one deduces from the above 
equation that they differ in our model, even without fluctuations of the scalar flelds: 




fHM'a 

Cv 



1 



z — 2" '^^iQ '^'^iQ ■ 

l-c^J Cv 



Appendix C: Calculation of the third-order action 



In this section, we give details about the calculation of the induced gravity action at third- 
order in the scalar perturbations, as needed in section (4.1). We work at leading order 
in the slow-varying approximation and neglect the metric perturbations (although it is 
straightfoward to include them), which is consistent in the regime k ~ 1. The expressions 
of the various quantities needed when starting from the ADM expression of the induced 
gravity action (in the second line of Eq. (3.6)) are as follows (the fleld space indices are 
contracted with the trivial metric Gjj = 6jj): 



(2) 



hij = a^6ij + fdiQ^djQi = a^6ij + 



(Cl) 
(C.2) 
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with 



With the notations 



= f^Q. (C.3) 
/3f ^ = fQAQ' (C.4) 



a« = (C.5) 



(2) 



/ 



2r2 



Q/Q' (C.6) 



we also have 



2 



- ^(aW)^^ + 3a«a(^) + - ^^^^ + O(g^) (C.7) 

VJi = a^(^l + l^dQ'dQi^ + 0{Q') (C.8) 

f;, = I^d'Q^d.dkQi + 0{Q') (C.9) 

^(3) = I_ {d'^Q'd^Qi - d.djQ'd'd^Qj) + O(g^) (C.IO) 

4- = a'H6,, - + - di^Sf^ + ^d.d.Qjd'^Q'dJ^'^ + 0{Q^) (C.ll) 

^ a 

f 

h'' = ^ - -.d'Q'd'Qi + 0{Q^) (C.12) 



and 



?(2) 



- ^^r^ (a^g.ag^ag/ + d.,d,Q^d'Q'd^Qi) + o(g^) . (c.13) 

We then have all the building blocks to calculate the second and third-order action 
in this regime. For the former, this straightforwardly gives the result Eq. (3.31). 
More work is required for the latter. In particular, to deal efficiently with the term 
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(d'd^^^^^ - 5*^(92/3(1)) (hlf - 2d(iP^f^ appearing in Eq. (C.13), it is useful to sepa- 
rate the scalar and vector contributions (in the sense of cosmological perturbations) from 
^l"^^ by writing 

^i'^ = f{d,S + W,) (C.14) 

with 

S = (dQ'dQi + Q'd^Qj^ (C.15) 

and d^Wi = 0. One then obtains 

_ J_ (^d'd^^^^) - 5^^-92^(1)) (/ig) - 2a(,^(f ) 

= ^ [o'^Q^dQ'dQi - d'd^Q^diQ'djQj + d'd^Q^d.djS - {d^Q^)d^s') .(C.16) 

Integrating by parts to make appear d'^S and using Eq. (F.2), one finds 

J dtd^x^ (^d^Q^dQ'dQj ~ d'd^Q^diQ'd.Qi + d'd^Q^didjS - {d^Q^)d^S^ 

= j dtd?x^- (^^dQld^Q^ + Qsd'd^QAdjQa-Qsd^Qsd^Qa^ • (C.17) 

Together with a few non-trivial integrations by part listed in appendix (F), and especially 
Eq. (F.5), one then obtains the result Eq. (4.2). 

Appendix D: Focusing on the induced gravity 

We have seen in our discussion of the background evolution in section 2.2 that the 
induced gravity should be subdominant compared to standard gravity - in the sense 
that ^ Cj)Mp and hence k ^ 1 - in order to achieve a phase of quasi de-Sitter 
inflation in the relativistic regime c|, ^ 1. Despite this, it is interesting to understand 
the opposite limit in which the induced gravity is dominant k ~ c|,. In this appendix, 
we therefore consider the model defined by the action (2.16) in the formal limit in which 
Mp — )■ 0, thereby focusing on the effect of the induced gravity only. We do not make 
a comprehensive study of this model bur rather try to make the links with some of the 
general results in the main part of this paper. 

In the main text, we argued that the single field second-order scalar action in terms 
of ( Eq. (3.16) can be recovered by going to the Einstein frame in which ^^j, is treated 
as a cosmological metric and applying the standard results of /c-inflation [20]. Indeed, in 
the uniform inflaton gauge - whose condition 6(j) = reads the same in both frames - the 
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spatial parts of the metrics g and g coincide (see Eq. (3.3)). Hence, the scalar degree 
of freedom in the Einstein frame C, defined such that hij = a^(t)e^^, actually coincides 
with (, such that hij = a^(t)e^''. As for the lapse functions A^ and A^ and the shift 
vectors Aj and Aj, they appear in both frames as auxiliary variables, and extremizing the 
action with one set or another leads to the same constraints. Therefore, for the purpose 
of determining the scalar action in terms of ( at any perturbative order, treating g or 
g as the metric of reference represents a mere change of variables in the intermediate 
calculations. 

The situation is more subtle in the multifield situation. Indeed, the standard variables 
to use then are the field fluctuations in the spatially fiat gauge. However, the two 
spatial metrics being related by hij = hij + fGijdi(l)^dj(f)'^ (Eq. (3.3)), the spatially fiat 
gauge for the cosmological metric does not coincide with the spatially fiat gauge for the 
induced metric, because of the spatial gradients of the scalar fields. However, the two 
gauges coincide at linear order in the fluctuations. Hence, expressing the brane action as 
'S'brane = / d!^ X F {X^^ , (j)) where = —^g'^'^d^cp^d^cf)'^, one can readily apply the 
linear analysis of Ref. [35] for such scalar field Lagrangians without making any change 
of variables. One then easily recovers the Mp — )■ limit of the result Eq. (3.21). In 
particular, note that the speed of propagations of scalar perturbations with respect to 
cosmic time is then unity. 

From the analysis of the linear cosmological perturbations, one deduces that the latter 
are not ghosts if and only if e > or equivalently, with 

-_3/a^ 1-fV 

2 CD 1 - CD + CvfV ' 

1 - fV > ^ 1 - 3a > 0. Given that Eq. (D.l) is equivalent to 

a + 3Ha + ^-^V,^ = 0, (D.2) 

where = ^iVi, one sees that the fields tend to climb up the potential precisely in the 
regime in which the fluctuations are ghosts. Note also that the deceleration parameter e 
is given by e — s, so that one can imagine, if s = cx)/{Hcx>) is sufficiently large, a ghost- 
free regime (e > 0) violating the null energy condition e < and in which the speed of 
propagation of all scalar fluctuations is unity (this should be contrasted with the results 
in Ref. [152] for example). Finally, note that if s can be neglected, e ~ e can be small 
in the ghost-free relativistic regime only if 1 — fV ^ c©, in other words, at the expense 
of having very small kinetic terms for the scalar fields and hence at the risk of being in a 
strongly coupled regime. 



(D.l) 
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Appendix E: Effects of the warping 

In this section we discuss some of the consequences of considering a general non-constant 
warp factor h in our starting point action Eq. (2.4). Using the general gravitational (2.10) 
and fields (2.14) equations of motion, and employing the same computational trick as in 
section 2.2, we easily determine the modified Friedmann equations 



and 



M2 



Pbr 



v + 



a2 



2c 



D 



M2 

CD 



f 



1 

CD 



1 



-D 



hi 



hi 



Cv 



H 



(E.l) 



(E.2) 



where H, defined by 



H = H 



4/ 



(E.3) 



differs from the Hubble parameter due to the non-trivial warping. As for the equations of 
motion for the scalar fields (2.14), they read, in the background: 



+ 



3if' 



M2 



,M2 



cvf 



{cUei+ 




)J,J ^ pCD 



2H + 3H^ + 2 



(E.4) 



where, like in the main text, = /(^ is the unit vector (with respect to the field space 
metric) pointing along the background trajectory in field space, we have introduced the 
projector on the subspace orthogonal to it 



-ij 



G 



IJ — CcrlCcrJ 



(E.5) 



and 



p^0/^0/ + r^^^0^0^ (E.6) 

represents the acceleration vector of the fields in curved coordinates (in field space) [153, 
154]. Because of the non-trivial warping, second-order derivatives of the fields appear not 

only through cv in Eq. (E.2) and through Vt<p^ in Eq. (E.4) but also in H, which contains 



f = f,i^' + lij^ 



(E.7) 
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Hence, should one wish to express separately and explicitly the second-order derivatives 
of the scale factor and of the fields, one would first need to solve a system of three linear 
equations for the three unknows H, a = e^jVtv ^"^^ -L/ f,j1^tv'- Eq- (E.2) and the 
contractions of Eq. (E.4) with e^j and _L/ f^j. This can easily be carried out but the 
explicit solutions are not particularly illuminating and we do not reproduce them. On 
the other hand, in the quasi de-Sitter slow-varying relativistic regime described in section 
(2.2), the equations of motion for the fields take the simple approximate form: 

3Ha{l-3a) + CT,V^-0, (E.8) 

where the induced gravity effects appear in the same combination 1 — 3a as in the grav- 
itational equation (2.28). Eventually, let us note that the same complexity as described 
above holds of course at the level of the perturbations: due the non-zero gradients of 
the warp factor along the entropic directions in field space, second-order time derivatives 
of the adiabatic and entropic perturbations do not decouple in the equations of motion, 
contrary to Eqs. (3.27) and (3.28). Hence, although there is no conceptual obstacles to 
deriving the exact second-order action in that case, the result would be of little practical 
analytical use. 



Appendix F: Useful integrations by part 

We collect here a number of non completely trivial integrations by part needed in the 
calculation of the third-order action in appendix C. 



dt d^x - ( Q^didjQ^d'd'Q^ - Q.{d'Q,] 



&t d^x - ( d'Q^dQ^dQ^ - d,Q^d.jQ^d'&Q^ 



1 fl 



a \2 



(F.3) 

j dtd?x ^ d,djQ,d'Qsd'Qs = j dtd'x ^ (^-d^QsdQ^dQs + ^d^QAdQsf^ (FA) 
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dtd'x - 2Q,d'd'Qsdid,Q, - 2Qsd'Qsd'Q„ + -{dQsYd'Q, + dQ.dQsd'Q 



(F.5) 
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